Quasi-elastic neutron scattering in the high-field phase of a Haldane antiferromagnet. 
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Inelastic neutron scattering experiments on the Haldane-gap quantum antiferromagnet 
Ni(C5Di4N2)2N3(PFe) are performed in magnetic fields below and above the critical field H c at 
which the gap closes. Quasi-elastic neutron scattering is found for H > H c indicating topological 
excitations in the high field phase. 



As first realized by Haldane, integer-spin one- 
dimensional (ID) Heisenberg antiferromagnets (AFs) are 
exotic "quantum spin liquids" with only short-range spin 
correlations and a gap in the magnetic excitation spec- 
trum [[l). Haldane-gap systems exhibit numerous un- 
usual properties, and particularly interesting predictions 
were made for their behavior in high magnetic fields 
@j u i Q i 0- ^ n an external field the gap excita- 
tions, which are a S = 1 triplet, are subject to Zeeman 
splitting. The gap A for one of the three branches de- 
creases with field (Fig. |jl(a)) and closes at some critical 
field H c ~ A/gnB @, |, |, !■ Near H c this problem 
is equivalent to Bose condensation || |j. However, at 
H > H c the formation of a condensate is prevented by 
strong interactions between magnons. As a result, the 
high-field phase is rather unusual, with power-law spin 
correlations, and is an example of a Luttiner spin liquid 
H M . The basic physics at H > H c can be understood 
within a simple model. The magnon interactions are as- 
sumed to be a hard-core repulsion, in which case the ex- 
citations behave as free fermions || |[ ||, |l(| . At H > H c 
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FIG. 1: Dispersion in the lowest-energy Haldane-gap mode 
below (a) and above (b) the critical field H c (blue). In the 
free-fermion model, excitations with a negative energy at H > 
H c form a Fermi sea (cyan) with a characteristic Fermi vector 
qf Lowest energy excited states are viewed as "particles" for 
|<7U — 7r| > <jf or "holes" for |<jy — 7r| < <jf (red line). 



the ground state is a Fermi sea of excitations with field- 
dependent Fermi density and wave vector qp, as illus- 
trated in Fig. 0(b). The spectrum is gapless, and dy- 
namic spin correlations are incommensurate. The effect 
is related to field-induced incommensurability in gapless 
spin chains, that is described using a different fermion 
mapping |ll| . 

For a long time the spin dynamics of Haldane gap AFs 
in the high-field phase evaded direct experimental inves- 
tigation by neutron scattering techniques. This was pri- 
marily due to large values of H c for most known model 
quasi-lD integer-spin compounds. The situation changed 
with the discovery of Ni(C 5 Di4N 2 ) 2 N 3 (PF 6 ) (NDMAP), 
a material well-suited for neutron scattering experiments 
and a critical field of only 5 T|l^]. For H > H c resid- 
ual interactions between the S = 1 Ni 2+ -chains in this 
compound lead to commensurate long-range AF order- 
ing at low temperatures u2, Ol. However, the strength 
of effective inter-chain coupling is about 0.2 meV, small 
compared to the in-chain exchange constant J « 2.8 meV 
fil] , p^| . At fku > 0.2 meV one can thus directly study the 
exotic dynamics of the ID model. In the present work we 
report inelastic neutron scattering studies of low-energy 
excitations in NDMAP, for magnetic field below and well 
above H c , at temperatures above the 3D-ordered phase. 

Experiments were carried out on the SPINS 3-axis 
spectrometer installed at the cold neutron facility at the 
National Institute of Standards and Technology Center 
for Neutron Research. 23 deuterated NDMAP single 
crystals were co-aligned to produce a sample of total 
mass 1.2 g and a mosaic of about 7°. The sample was 
mounted in a cryomagnet and measurements were per- 
formed at T = 2.5 K and magnetic fields H = 0, 2, 4, 5.5, 
7 and 9 T. The field was applied along the b axis of the 
orthorhombic structure (space group Pnrnri), perpendic- 
ular to the spin chains that run along c (cell parameters 
a = 18.05 A, b = 8.71 A, c = 6.14 A). The spectrometer 
employed a vertically focusing pyrolitic graphite (PG) 
(002) monochromator and a 23 cmxl5 cm flat multi- 
blade PG(002) analyzer positioned 91.4 cm from the sam- 
ple. Scattered neutrons were registered by a position- 
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FIG. 2: False color plots of inelastic intensity measured in 
NDMAP at T = 2.5 K for different values of magnetic field 
applied along the b axis. The extra intensity seen below 
q« = 0.4c* at low energies in (c) is most likely an artifact 
of imperfect background subtraction. Solid lines are as de- 
scribed in the text. 



sensitive detector (PSD). Measurements were done with 
fixed-final neutron energy for the central analyzer blade 
of 3.1 meV. The collimation setup was ( 58 Ni-guide)- 
(open)-(open)-80'(radial), and a Be filter was used ei- 
ther before or after the sample. To minimize the contri- 
butions of the sample mosaic to wave vector resolution, 
all scans were performed such that momentum transfer 
for the center pixel of the PSD was along c* . Wave vector 
transfer along c* will be denoted by qii . In a single set- 
ting of the spectrometer the PSD covered an energy range 
of 1 meV and a range of wave vector transfer along the 
chain of 0.3 c* Separate scans were corrected for the ana- 
lyzer efficiency and merged to produce a single data set. 
A linear background was subtracted from each const-S 
slice taken from the data set at each field. The data were 
smoothed using a Gaussian kernel with SE = 0.05 meV 
FWHM and <5<7|| = 0.01 c*. Typical resulting data sets 
are visualized in Fig. 0. 

For the following discussion it is important to estimate 
the critical field H c . Long-range ordering, that is es- 
sentially a three-dimensional (3D) effect and occurs at 
Hl D « 12 T for T = 2.5 K, is of little interest in the 
present study. The relevant quantity is the field H c at 
which the gap for the lower-energy Haldane excitation 
vanishes at = (0, 0, 0.5), where all the measurements were 
performed. H c can be estimated from the gap energies 
measured in zero field |14j and the known gyromagnetic 
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FIG. 3: Field dependence of the gap energies (a) and ex- 
citation intensities (b) measured in NDMAP at T = 2.5 K 
(symbols), (c) A constant-Q scan collected at H = 5.5 T 
(symbols). Solid lines and shaded areas are as described in 
the text. 



ratio for the S = 1-carrying Ni 2+ -ions in NDMAP [J12J. 
Using Eq. 2.14 in Ref.@ we get H c w 7.3 T. Below we 
shall first review the results obtained for H < H c . and 
proceed to discuss the 9 T data, that represent the high- 
field phase and contain our most important findings. 

The Haldane triplet in NDMAP is split by single-ion 
anisotropy [Q. The parabolic shape seen in Fig. ||(a) 
for H — corresponds to two lower-energy gap excita- 
tions polarized along the a and b axes, with gap ener- 
gies Ai = 0.48(1) meV and A2 = 0.65(1) meV, respec- 
tively. The experimental resolution is not sufficient to 
resolve these two modes. The gap in the c-axis polar- 
ized Haldane excitation, A3 ps 1.9 meV [14j, is outside 
the energy range covered in the present experiment. In 
Fig. 0(b), at H = 5.5 T, one clearly sees the splitting of 
the doublet. Excitations polarized along the b axis (up- 
per branch) are S z = states, z being chosen along the 
applied field. To first order, they are not affected by the 
magnetic field. The gap for the a-axis polarized Haldane 
excitation (lower mode) decreases with field by virtue of 
the Zeeman effect. 

The data were analyzed using a parameterized model 
cross section, numerically convoluted with the calculated 
spectrometer resolution. Details of the data analysis will 
be given elsewhere and only the main points are sum- 
marized here. The cross section was written in the sin- 
gle mode approximation fl4J ]. For each observed ex- 
citation branch the dispersion relation was written as 
(fiujq) 2 = A 2 +v 2 sin 2 (qc) (Model 1). The gaps and inten- 
sities for each branch were adjusted to best fit the experi- 
mental data. At each field the entire data set with energy 
transfers between 0.2 meV and 1.6 meV was analyzed 
with a single set of parameters in a global fit. The spin 
wave velocity v was fixed at v = 6.9 meV, as determined 



at H = 0. The model was found to describe the data 
measured at up to 7 T very well, with a residual % 2 < 1.5 
in all cases. The obtained field dependence of the gap 
energies is plotted in Fig. ||(a). The corresponding dis- 
persion relations are shown in solid lines in Fig. 0(a)-(b). 
Fig. f|(c) shows a constant-Q scan extracted from the 
5.5 T data set (symbols) plotted over the profile simu- 
lated using the refined parameter values (black line). The 
shaded areas show partial contributions from the two gap 
modes. Similarly, Fig.H(a)-(c) shows a series of extracted 
constant-I? scans for H = 5.5 T, plotted on top of the cor- 
responding simulated profiles (red lines). The measured 
field dependence of gaps follows the general trend [solid 
lines in Fig. |3|(a)] expected from a perturbation treat- 
ment of the Zeeman term |l6|]. In agreement with the 
perturbation analysis [blue line in Fig. |3|(b)], the mea- 
sured intensity of the upper mode is field independent 
[Fig. ||(b), cyan symbols] while that of the lower mode 
decrease with H as it acquires polarization along the c- 
direction to which our experiment is insensitive [Fig. |3j(c), 
green line]. At H < H c the behavior is well described 
by first order perturbation theory and is similar to that 
previously found in the extensively studied Haldane gap 
compound Ni(C 2 H 8 N 2 )2N02(C104) (NENP) 0. 

We now turn to the data visualized in Fig. ||(c), col- 
lected at H — 9 T> H c . The gap in the upper mode 
has increased substantially, to about 0.93 meV, as com- 
pared to 0.63 meV at H = 5.5 T. Interestingly, below this 
gap, a weak but well-defined area of scattering extends 
all the way down to the lowest accessible energy trans- 
fer. Due to limited wave vector resolution and the steep 
1/lo intensity scaling, the dispersion of these excitations 
is difficult to discern other than by a quantitative anal- 
ysis (see below). This is not conventional critical scat- 
tering associated with the imminent phase transition, as 
the experiment was done away from Q c = (h, 0.5, 0.5) 
where magnetic Bragg peaks appear below Tjv- This 
first experimental observation of quasi-one-dimensional 
quasi-elastic neutron scattering for H > H c implies that 
individual spin chains in NDMAP are in a critical phase 
from which quasi-long-range order would likely emerge 
on cooling even in the absence of inter-chain coupling. 

Closer examination of the data reveals important 
changes at high fields beyond closing the gap in the neu- 
tron scattering spectrum. Applying the low field analysis 
procedure to the 9 T data yields a significantly worse fit 
(x 2 = 2.0) than was achieved for H < H c . The prob- 
lem is mainly in the wave vector and energy dependence 
of inelastic scattering at low energies. The discrepancy 
is particularly pronounced in the shape of constant-.^ 
scans, where the data shows peaks considerably broader 
than simulations based on parameters obtained in the fit. 
To quantify this behavior we scaled the measured inten- 
sity by the energy transfer and integrated the result over 
the range 0.25-0.5 meV. The resulting g-dependence was 
fit to a Gaussian profile to obtain an intrinsic (resolution 
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FIG. 4: Symbols: constant-i? scans collected in NDMAP at 
H = 5.5 T< H c (a-c) and H = 9 T> H c (d-f). Lines are as 
described in the text. Arrows indicate the projection of the 
resolution ellipse onto the chain axis. 



corrected) width of 0.048(5)c*. The peak simulated using 
parameters determined in the global fit has an intrinsic 
width of only 0.023c*. The anomalous q- width is also 
visible in each measured constant- E scan individually, as 
for example at 0.4 meV energy transfer [Fig. |j (d)]. It 
is important to emphasize that the discrepancy is not 
related to an error in resolution calculation. Indeed, at 
lower fields, e. g. at H — 5.5 T [Fig. [| (a)] Gaussian 
fits to const-i? scans and simulated profiles yield identi- 
cal widths. Moreover, even at H — 9 T at higher energy 
transfers, where the dominant contribution is from the 
0.93 meV mode, the model cross section fits the data 
quite well [Fig. |] (e) and (f), solid lines]. 

At least two theoretical frame works could potentially 
account for the gii —broadened quasi-elastic scattering in 
the high field phase of NDMAP. A classical easy plane 
antifcrromagnet with a field in the easy plane can be 
mapped onto a Sine-Gordon model with a gapped spin 
wave spectrum and topological soliton excitations. |l8| 
While solitons have a finite rest mass, they are sta- 
ble topological objects that move freely along the spin 
chain. The soliton gas yields quasi-elastic neutron scat- 
tering with a <7|| —width that varies in proportion to 
the soliton density. The soliton model was used suc- 
cessfully to account for the phase diagram of NDMAP 
with fields normal to the chain axis.jjli In addition Soli- 
ton pinning at defects may account for the quasi-two- 
dimensional frozen state in NDMAP for fields perpen- 
dicular to the easy plane. Q Further evaluation of the 
Sine-Gordon model for the high field phase of NDMAP 
would require measurements of the H/T dependence of 



the g||— width for comparison to the theoretical soliton 
density. The free-fermion model presented in the intro- 
duction describes a different type of topological excita- 
tions. In this model at H > H c the lowest energy exci- 
tations with momentum transfers along the chains near 
<7|| = ir/c are single-particle or single-hole states. These 
have a dispersion relation as shown by the red line in 
[Fig. p] (b)]. To compare the 9 T data to the free fermion 
model, we constructed an empirical dispersion relation 
for the lower mode [Model 2] with qp as a new param- 
eter. The dispersion was assumed to be linear outside 
the region ir/c — qp < q» < n/c + q-p, and to follow half 
a period of a sine wave within these boundaries. The 
prefactor for the sinusoidal part was chosen to give the 
correct spin wave velocity at gii = tt/c ± qp. As illus- 
trated in Fig. [I] (b), a given value of qp corresponds to 
a "negative gap" A = -2q F cv/ir. With q F = 0.03(l)c* 
(A = —0.3(1) meV) Model 2 yields a considerably bet- 
ter agreement with experiment (x 2 = 1.3) than Model 
1. The refined gap energy for the upper mode and the 
"negative gap" for the lower branch are indicated with 
triangles in Fig. (a) . Simulations based on resulting pa- 
rameter values are shown in magenta lines Fig. (d)-(f), 
and the dispersion relations are shown in a solid line in 
Fig. ||(c) . While at high energies where the upper mode 
is dominant the two Models are virtually indistinguish- 
able, below 0.6 meV transfer Model 2 does a better job 
at reproducing the observed peak widths [Fig. || (d)]. 

The free fermion model makes a very robust predic- 
tion for the field dependence of qp- The Fermi vector 
q F determines the number of S z = 1-carrying particles 
present in the ground state at H > H c and is there- 
fore directly proportional to the uniform magnetization: 
M = g/iB 2q F /c* , where M is the magnetization per spin. 
For NDMAP the magnetization curve has been mea- 
sured experimentally [g0| and at 9 T M ~ 0.1/XB/Ni 2+ . 
This gives qp = 0.025c*, remarkably close to the exper- 
imental value. It is important to stress that the free- 
fermion model is a rather crude approximation. In fact, 
at H > H c the system is to be described as a sea of 
interacting fermions, i. e. a Luttinger liquid (for a re- 
cent discussion see Ref. 0). As a result of these inter- 
actions, there is no well-defined Fermi surface or even 
single-particle poles in the dynamic susceptibility. In- 
commensurate correlations are obscured: the equal-time 
spin correlation function is a broad peak at the com- 
mensurate q\\ = tt point. It is probably for this reason 
that static long-range ordering in NDMAP [O and the re- 
lated compound Ni(C 5 Hi4N2)2N 3 (C104) (NDMAZ) @ 
is commensurate. In the Luttinger liquid model the low 
energy excitations are a continuum. Experimental evi- 
dence for this behavior was recently obtained in studies 
of NDMAP [0 and NENP (22). While the resolution of 
our experiments on NDMAP is insufficient to distinguish 



for the data analysis. It can reproduce the measured data 
quite well, and yields a self-consistent estimate for qp. 

In summary, our experiments demonstrate 
q» —broadened quasi-elastic scattering in the high- 
field phase of a S = 1 Haldane-gap compound. Various 
models based on topological excitations with low energy 
phase fluctuations may account for the data, including 
the Sine-Gordon and the free fermion models. Further 
theoretical and experimental work will be required to 
identify the correct description of the high field phase in 
the anisotropic spin-1 chain. 
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